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Abstract 

We discuss the tension between discrete flavour symmetries and extended scalar 
sectors arising from lepton flavour violation experiments. The key point is that 
extended scalar sectors will generically lead to flavour changing neutral currents, 
which are strongly constrained by experiments. Because of the large parameter space 
in the scalar sector such models will, however, usually have no big problems with 
existing and future bounds (even though the models might be constrained). This 
changes considerably once a flavour symmetry is imposed in addition: Because of the 
symmetry, additional relations between the different couplings arise and cancellations 
become impossible in certain cases. The experimental bounds will then constrain the 
model severely and can easily exclude it. We consider two examples which show how 
these considerations are realized. The same logic should apply to a much wider class 
of models. 
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1 Introduction 



The Standard Model (SM) is a very successful theory. Apart from only describing the 
phenomena observed, it has even made predictions of new particles such as the Z-boson or 
the t-quark. There is only one missing piece, namely the well-known neutral scalar Higgs 
boson, that will hopefully be found at the LHC collider. 

Apart from this success there are, however, some observations which indicate that the SM 
is incomplete, among them, e.g., the observation of Dark Matter [T], the baryon asymme- 
try [2] , and the hierarchy problem [3] . An extension of the SM is therefore necessary (this 
is usually called Physics beyond the SM (BSM)). The various possible extensions of the SM 
very often contain extended scalar sectors. Other extensions address the flavour problem 
and introduce new flavour symmetries to explain the apparent regularities of masses and 
mixings. 

In this paper we study the difficulties arising when one tries to combine a model with an 
extended scalar sector with a discrete flavour symmetry. The key point is that there are 
actually quite strong constraints on models with extended scalar sectors. Since, however, 
the scalar sector of a theory is in most of the cases poorly known (meaning that there are a 
lot of free parameters), such a model can usually not be excluded easily, because of internal 
cancellations between several of the parameters that may cause some observables to nearly 
vanish. If, on the other hand, some additional structure is imposed on the model (by, e.g., 
a discrete flavour symmetry), then additional relations between some of the parameters 
can easily rule out the corresponding model or to at least restrict its parameters to very 
narrow ranges. 

The paper is organized as follows: In Sec. [2, we introduce the argumentation which lead 
us to our statement that models with extended scalar sectors may get into trouble by the 
introduction of an additional flavour symmetry. This is exemplified in Sec. [31 where we 
present two particular models for which our logic clearly works. The numerical results that 
we have obtained are presented and discussed in Sec. IH and we finally conclude in Sec. [51 
The basic properties of the discrete groups that we use are given in the Appendices A (A±) 
and B (D 4 ). 

2 The general arguments 

A natural way to extend the SM is to add further scalar particles, which have not yet 
been discovered. These could, e.g., be additional 5't/(2)-singlets [I], doublets ("Two Higgs 
doublet model", THDM), or triplets j5]. Depending on the model, it can then be the case 
that more than one Higgs field contribute to the masses of all particles or that certain 
Higgses only give masses to a particular choice of particles [6]. These models will then, 
however, generically lead to flavour changing neutral currents (FCNCs) [7j and hence to 
lepton flavour violation (LFV) processes [8], which are quite strongly constrained [9]. It is, 
however, also not easy to rule them out that way, since they will in general yield complex 
3x3 Yukawa coupling matrices, which hold a lot of freedom in their 18 parameters. So, in 
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most of the cases, such a model will be able to fit all neutrino data without any problems, 
even if it is strongly constrained. 

On the other hand, there are also ways to impose more structure onto the SM in order 
to get an understanding of quantities like mixing angles, or so. This is usually done by 
so-called (discrete) "flavour symmetries" under which the SM-fermions (and, depending on 
the model, also (additional) scalars) are charged in a certain way. If, e.g., two generations 
of SU{2) L doublets are components of the same doublet representation of a discrete flavour 
symmetry (such as the dihedral groups D 3 ~ S3 [10] or D 4 [HI E2l [13]), then this property 
will generically lead to /x-r symmetry [H] by which two mixing angles are predicted, #23 = 
7r/4 and #13 = 0. Moreover, in order to increase the predictivity, one can also assign all 
three generations of SU (2)^ doublets to form a triplet of a discrete flavour symmetry (such 
as A 4 [15j ESI [T7J EH EH]). This can lead to tri-bimaximal mixing [20] . in which also #12 is 
fixed to be tan 9 12 = 1/ \/2. 

Imposing such symmetries adds more structure to the model in the sense that one obtains 
relations between different entries of the Yukawa matrices. By that way, one can obtain 
the neutrino oscillation parameters as well as the charged lepton masses as functions of 
only a few parameters, which can then be checked on whether they are in accordance 
with data, or not. However, such models generically need a lot of scalars in order to 
break the flavour symmetry in a valid way. In case the normal Higgses are not charged 
under the flavour symmetry, these are additional SM-singlet scalars ("flavons"), which are 
only charged under the discrete symmetry and can hence break it by obtaining a vacuum 
expectation value (VEV). These scalars will, again, lead to horribly large FCNCs, which 
crashes with phenomenology. 

One way out is to decouple the flavons by giving them masses associated with the breaking 
scale of the flavour symmetry, which can be much higher then the electroweak scale. This 
is, of course, somehow only hiding the problem, but it will make the model fit better. 

We now apply the following logic: 

1. We impose a flavour symmetry and decouple the flavons in order to end up with an 
effective low energy model with a scalar sector that is slightly extended compared to 
the SM. This could, e.g., be a THDM or something similar. 

2. This procedure should make the model fit better, because the possible problems that 
could arise by the flavons are avoided. 

3. Since we have gained predictivity by imposing the flavour symmetry, we can fit the 
model to neutrino data, which allows us to extract certain ranges for the model 
parameters. 

4. The model, however, still has additional scalars compared to the SM, which will be 
able to mediate LFV-processes, whose branching ratios can be predicted using the 
fitted parameter values. 
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5. If this prediction does not fit with present (future) LFV-bounds, we are (will be) able 
to exclude the particular flavour symmetry imposed (in a certain scenario). Note 
that this logic will also hold in the non-decoupling case if no extreme fine-tuning is 
involved. 

In principle, this could work for any model with a slightly extended scalar sector. If the 
structure of the model is not extremely peculiar, which is rarely the case in the scalar sector 
of a theory, the additional scalars (compared to the SM) will unavoidably lead to LFV- 
processes, which are already strongly constrained. The key point is that these constraints 
are so strong, that imposing some more structure by adding a flavour symmetry can easily 
destroy the consistency of the model with all data. 

Here, we want to present such an analysis for one particular example, namely for Ma's 
scotogenic model [21], as this consists of a very minimal extension of the SM. Furthermore, 
it does not have too many possible LFV-diagrams, so that our logic is not shadowed by a 
heavy calculational apparatus. In this model, one can see immediately the effect of certain 
symmetries: Without imposing a flavour symmetry, one constrains quantities like 

\h* n h 21 + h* 12 h 22 + h* 13 h 23 \ (1) 

by LFV-processes like /i — > cy [22] , where h is the Yukawa coupling matrix involved. Such 
a combination can easily become zero for unfortunate values of some phases, exactly as 
the effective neutrino mass in neutrino-less double beta processes [23] . Imposing relations 
between certain elements of h hinders such cancellations to appear, and the term in Eq. (CQ) 
will generically be much larger than zero. 

We want to stress, however, that this particular model is just an example and that our 
idea may work for a much wider class of models. 

3 Constraining particular models 

3.1 One possible example: The scotogenic model 

There are a lot of different models for neutrino mass generation on the market [23]. A 
difficult task for all of them is to explain the smallness of neutrino masses compared to 
other particles we know in Nature. 

One way is to forbid a tree-level mass term for neutrinos and generate neutrino masses only 
by radiative corrections, as done in several models [211 [231 [251 [261 [271 [28]. Out of those, Ma's 
"scotogenic" model [21] (that we call "Ma-model" for simplicity) is particularly attractive: 
By adding only one additional Higgs doublet and heavy right-handed neutrinos to the SM, 
as well as imposing an additional ^-symmetry, it allows for sufficiently small neutrino 
masses. These masses are generated radiatively, because the additional neutral Higgs does 
not obtain a VEV that could lead to a tree-level neutrino mass term. Furthermore, due 
to the ^-symmetry, this model also provides a stable Dark Matter candidate, namely the 
lightest of the heavy neutrinos [29] or the lightest neutral scalar [30J. Constraints on the 
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model arise from various different sources as, e.g., lepton flavour violation or the Dark 
Matter abundance [22] • In that sense, this model is very "complete". 
The basic ingredients apart from the SM are: 

• 3 heavy right-handed (Majorana) neutrinos Nk, which are singlets under £77(2) and 
have no hypercharge 

• a second Higgs doublet r\ with SM-like quantum numbers that does not obtain a 
VEV 

• an additional ^-parity under which all SM-particles are even, while Nk as well as r\ 

are odd 

The corresponding Higgs potential looks like 

V = m?0V+m^ t 77+Y(0V) 2 +Y^ t ^)VA 3 (0V)(?7 t ?7)+A4(0 t 77)(77V)+Y Wvf + h.c] , 

(2) 

where <fi is the SM-Higgs. If m\ < and m 2 > 0, then only 0° will obtain a VEV 
v = 174 GeV, while (77°) = 0. Then, the Yukawa Lagrangian is given by 

£ Y = fijtf-Vi + (0°)%)e, c + hijtfvi - v + h)Nj + h.c, (3) 

which does not lead to a tree-level neutrino mass term due to the vanishing VEV of rj°. The 
neutrino masses can, however, be generated radiatively, which gives a natural suppression 
of the neutrino mass eigenvalues and can exploit the heaviness of the Nk (with masses Mk) 
as well. The mass matrix of the light neutrinos reads 

3 

(M u )ij = h ikhjk^k, (4) 
fc=i 

where 
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. f m 2 (H°) \ m 2 (A°) ( m 2 (A°) \ 



(5) 



_m 2 {H®)-M 2 K \ Ml J m 2 {A°)-M 2 K \ M 2 K J 

Note that we have named the Higgses like in the general THDM, with a = (3 = m\i = 
A 6>7 = [31J. The resulting Higgs masses are given by 

m 2 (h°) = 2\ lV 2 , m 2 (H°) =m 2 2 + (A 3 + A 4 + X 5 )v 2 , m 2 (A°) =m 2 2 + (A 3 + A 4 - A 5 )^ 2 , 
and m 2 (77 ± ) = m\ + X 3 v 2 . (6) 



3.2 The flavour symmetries considered 

In the following, we will present two models which constrain the structure of the Yukawa 
coupling matrix h in Eq. dSJ), without discussing a particular mechanism for vacuum align- 
ment. □ The first one, based on Ref. [32], represents the class of models which predicts 
tri-bimaximal mixing. The second one represents the class which predicts fi — r symmetry. 

■"■In general, the vacuum alignment can be achieved by a minimization of the scalar potential. 
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Table 1: The particle content of Model 1: The SM particles are 
the three left-handed lepton SU{2) L doublets l i} the right-handed 
charged leptons e?, and the SM-Higgs 0. The BSM particles are 
the right-handed neutrinos Ni, the second Higgs doublet t] (which 
does not obtain a VEV), and the flavons (ps, <fT, and Xi that only 
transform under v4 4 x Z^ nx . 



3.2.1 The A 4 -model (Model 1) 

The particle content of this model is given in Tab. [TJ The Lagrangian which is invariant 
under the flavour symmetry x Z| ux readg^l 



+ytj{ktfT2 + k<PT3 + hVTi)e c x + | [j/i[(2ZiiVi - / 2 iV 3 - hN 2 )^si 
+ (2l 3 N 3 - hN 2 - hNjipsa + (2^ - hN 3 - hNj^] 
+y 2 (l 1 N 1 + l 2 N 3 + l 3 N 2 ) x + M{N l N l + + N 3 N 2 ). 

Let us assume that the flavons obtain their VEVs as follows, 



(7) 




w s 






and (x) = u 



(8) 



and the SM Higgs gets the VEV (0) = v. Then, the Yukawa coupling matrix and the 
right-handed neutrino mass matrix for Model 1 can be written as 



'2a + b -a -a \ A (T 

-a 2a b - a and M R = M 1 
-a b- a 2a I \0 1 0, 



(9) 



where a = yi^r- and b = y 2 \- 

The charged lepton mass matrix in this model is diagonal, 



m M = -^y 2 w T , m 7 



A 



y 3 w T . 



(10) 



2 Here, we neglect the anti-symmetric part of the coupling between I and N or assume that the anti- 
symmetric coupling vanishes, which is done similarly in Refs. [321 133] . 
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Here, the hierarchies in the charged lepton masses are determined by the Yukawa couplings. 
Assuming that the Yukawa coupling of the r, x/f, is of 0(1) and the Higgs VEV v is 174 GeV, 
we can determine the ratio of the flavon over the cutoff scale A (^-) as being of the order 
of the Cabibbo angle squared, A 2 ~ 0.04. 

In order to make the discussion easier, we go to the basis where the right-handed neutrino 
mass matrix is diagonal. The matrix MrM r is diagonalized by the unitary matrix U r 



U r = 1 . (11) 




Note that the right-handed neutrino masses are degenerate, Mi j2)3 = M. 
The Yukawa coupling in this basis reads 

-a —a 2a + s 
\b-a 2a -a 
2a b — a —a 

Using Eq. (J4]), the neutrino mass matrix can be written as 

'(6a 2 + Aab + 6 2 ) -a(3a + 26) -a(3a + 26) 
M v = A lj2>3 \ -a(3a + 26) (6a 2 - 2a6 + 6 2 ) a(-3a + 46) (13) 
-a(3a + 26) a(-3a + 46) (6a 2 - 2a6 + 6 2 ) 

where Ai .2,3 = Ai = A 2 = A 3 , and the neutrino masses are given by the eigenvalues of 
M v Ml: 

m\ = (3a + 6) 4 A 2 j2i3 , m 2 = 6 4 A 2 >2j3 , and m 2 = (-3a + 6) 4 A 2 j2i3 , (14) 

which correspond to the eigenvectors (— 2, 1, l) T /v^6, (1, 1, l) T /\/3, and (0, — 1, l) T /\^2, 
respectively. In this model, the neutrino masses obey normal mass ordering. 
The neutrino mixing observables look like: 

Am| = (6 4 -(3a+6) 4 )A 2 2 3 , Am 2 A = -24a6(9a 2 +6 2 )A 2 2 3 , tan^ 2 = 9 13 = 0, and #23 = 

V2 

(15) 

In this model, we have only three free parameters (a, 6, M) to fit all observables. Therefore, 
this model is quite predictive (and hence harder to fit). 

3.2.2 The D 4 -model (Model 2) 

The particle content of this model is given in Tab. [2j The Lagrangian which is invariant 
under the flavour symmetry _D 4 x Z| ux reads 

Ci = ytheZjVe + yZikeZ + heQj^pe + yZfeeZ-he&jXe 

+y 1 l 1 N 1 ^ip v + y 2 Mi + hth)Nij + y 3 (h^2 - h^i)N 2 j + y^h^i - hifa)N s ^ 

+ ^M 1 N 1 N 1 + ^M 2 N 2 N 2 + ^M 3 N 3 N 3 . (16) 
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Table 2: The particle content of Model 2: Th SM particles are 
the three left-handed lepton SU(2)i doublets k, the right-handed 
charged leptons ef, and the SM-Higgs <fi. The BSM particles are 
the right-handed neutrinos JVj, second Higgs doublet 77 (which does 
not obtain a VEV),and the flavons <p e , Xe, <Pv, and V'i? that only 
transform under _D 4 x Z| ux . 



Let us assume that the flavons obtain their VEVs as follows: 

(ip e ) = u e , (xe) = -w e , (ip v ) = u , and ( j^j J = w ( _^ J , (17) 

and the SM Higgs gets the VEV (0) = Then, the Yukawa coupling matrix for Model 2 
can be written as 

'a N 
b — c d 
—b — c d> 

where a = y ± f, b = y 2 f, c = y 3 f, and d = y 4 f ■ 

The charged lepton and right-handed neutrino mass matrices in this model are diagonal, 

V V V 

m e = J^ytUe, «V = jivlUe - y^W e ), 7Tl T = -(y 2 U e + ytw e ). (19) 

Here, the hierarchy between the masses of e and (/x, r) arises from the smallness of the 
yukawa coupling yf. As we did for Model 1, we assume that the ratio (^-) is of order 
A 2 ~ 0.04. 

Using Eq. (J4]), the neutrino mass matrix can be written as 

(a 2 A 1 abAx -abA 1 

abA 1 b 2 A x + c 2 A 2 + d 2 A 3 -b 2 A x + c 2 A 2 + d 2 A 3 | . (20) 
-abAi -b 2 A 1 + c 2 A 2 + d 2 A 3 b 2 A 1 + c 2 A 2 + d 2 A 3 

The neutrino masses are given by the eigenvalues of M v Ml, 

m\ = 0, m\ = (a 2 + 26 2 ) 2 A 2 , and m\ = 4(c 2 A 2 + d 2 A 3 ) 2 , (21) 

which correspond to the eigenvectors 

{2b/ a, -1, 1) T , = h n = (-b/a, 1, 1) T , and (0, 1, 1) T /V2, (22) 
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Scenario 


m(h°) 


m(H°) 


m(A°) 




a 


120.0 


32.9 


84.5 


93.0 


(3 


120.0 


60.4 


101.5 


111.5 


7 


120.0 


946.8 


950.0 


950.3 


5 


120.0 


548.9 


549.4 


550.6 



Table 3: The Higgs masses (in GeV) for the different scenarios defined in Eq. ( 127]) . 



Quantity 




{Am\) aoT , 


#12 


#13 


#23 


Best-fit 


7.67- 10~ 5 eV^ 


2.46 • 10~ 3 eV^ 


34.5° 


0.0° 


42.3° 


la 


2.15 • 10~ 6 eV^ 


0.15 ■ 10" a eV^ 


1.4° 


7.9° 


4.2° 



Table 4: The neutrino mixing parameters (best-fit values and symmetrized la-ranges) 
obtained by a global fit [54"] . 

respectively. 

In this model, the neutrino masses will obey normal ordering. 
The neutrino mixing observables look like: 

Am 2 & = (« 2 + 26 2 ) 2 A?, Am\ = 4(c 2 A 2 + d 2 A 3 ) 2 , tan# 12 = 6 13 = 0, and 9 23 = -. 

\/1b 4 

(23) 

In this model, we have 7 free parameters (a, b, c, d, Mi, M 2 , M 3 ) to fit all neutrino observ- 
ables. This makes Model 2 much easier to fit, but we of course pay the price of losing 
predictivity. 

3.3 Phenomenological analysis 
3.3.1 The general procedure 

In this section, we describe the analysis procedure we have applied. The first thing to say 
is that there are constraints that are required for a THDM like in Eq. ([2]) (Ai > 0, A 2 > 0, 
A 3 > — VA1A2, and A3 + A4 — | A5 1 > — \f\1\2] they keep the potential stable) as well as 
consistency conditions for a Ma-like model (m\ < and m\ > 0; these are necessary in order 
for 0° to obtain a VEV, while rj° obtains none). Furthermore, there are limits from direct 
searches at collider experiments [35J: m(h°) > 112.9 GeV and m(i? ± ) > 78.6 GeV, both at 
95% confidence levelo Further constraints arise from the W- and Z-boson decay widths, 
namely m(iy ± ) + m(H°), mtf*) + m(A°) > M w and 2m(iy ± ), m(H°) + m(A°) > M z , 
as well as from the requirement of perturbativity for the Higgs potential, A 2 < 1 and 
X 2 3 + (A 3 + A 4 ) 2 + \\ < \2X\ [30]. 

Strong constraints also come from the correction to the p-parameter [36]. The explicit 

3 Note that these constraints do not apply to the "inert" Higgses H° and A . They are constrained 
much less severely by the current limits, differently from a normal THDM. 



s 



formula for this correction reads 

A P = tP^JT ■ [Hm 2 2 ,m 2 (H )) + F(m 2 2 ,m 2 (A )) - F(m 2 (H°), m 2 (A ))] , (24) 
where 

C x±y xy_ i x t i y 

v yj 1 0, for x = y, v 1 



and ct(Mz) = 1/127.9. The experimental constraint is 

Ap = -0.0006 ±0.0008, (26) 

which cuts the allowed parameter space for the Ma-model. Since we want to focus on 
neutrino physics and lepton flavour violation, we do not try to fit the Higgs sector as well, 
but rather use four different benchmark scenarios that all fulfill the consistency conditions, 
as well as the experimental bounds from direct searches and from the measurement of 
the correction to the p-parameter (at 3a). In the form (mi, rri2, Xi, A2, A3, A4, A5), these 
scenarios are: 



a 

7 
5 



(lOOzGeV, 75GeV, 0.24, 0.10, 0.10, -0.15, -0.10) 

(lOOiGeV, 98.5GeV, 0.24, 0.30, 0.09, -0.18, -0.11) 

(lOOiGeV, 950GeV, 0.24, 0.50, 0.02, -0.12, -0.10) 

(lOO^GeV, 550GeV, 0.24, 0.30, 0.02, -0.05, -0.01) (27) 



The corresponding Higgs masses are given in Tab. [31 We have chosen these four scenarios 
such that they are consistent with the 3<r-range of WMAP-data for H° being the Dark 
Matter candidate, which cuts the allowed parameter space significantly [30]. This leads to 
some more consistency conditions, as H° has to be the lightest of all scalars and it also 
has to be lighter than the heavy right-handed neutrinos. 
For all these scenarios, we do the following: 

1. First, the models are fitted to neutrino oscillation data, i.e., mixing angles and mass 
square differences [HI]. This is done by the x 2 -function 



(ft - € p " 2 



N 

X 2 = E l " ^ . (28) 
i=i Gi 



where g« are the observables obtained from neutrino oscillations (# 12 , # 13 , 6*23, Am^, 
AttIq), which are calculated in terms of the model parameters (cf. Sec. 13.21) . g° xp are 
their measured counterparts and <jj are the corresponding (symmetrized) standard 
deviations. The best-fit model parameters are determined by a minimization of the 
X 2 -function. By projection on the different directions in the parameter space, we 
determine the la- and 3<r-ranges of the model parameters. 
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2. Next, we calculate the maximum and minimum values of the quantities measured in 
different LFV-experiments (/x — > e^y, r — > /X7, t — > e^y, and /x-e conversion for four 
different nuclei) by varying the model parameters within their la- and 3<r-ranges. 

3. Finally, we compare how well different past and future LFV-experiments are able to 
constrain or exclude the particular model in the four scenarios. 



3.3.2 The x 2 -fit 

After outlining the general points, we will explain the procedure in more detail using 
scenario a (cf. Eq. ( 1271) ) in connection with Model 1 (cf. Sec. 13.2. ip as example. 
The x 2 -function has already been given in Eq. fl28l) and the experimental values and errors 
of the neutrino observables are summarized in Tab. HI These observables in terms of model 
parameters have been given in Eq. (|15p . The minimization of the x 2 -function then yields 
the following best-fit values for the three parameters: 

a = 0.0189, b = -0.691, M = 2.42 • 10 6 GeV. (29) 



Note that the parameter b is negative to fit the normal mass ordering, cf. Eq. ( fl~5i) . In the 
minimization we have required M\^,z > fn(H°) and Mi^,z > Mz/2 for consistency reasons. 
The 1<T-(3(T-) values for the model parameters are obtained by inserting all values from 
Eq. (|29|) into the x 2 -function, except for the one parameter that is to be constrained, and 



by determining the intersections of the remaining 1-dimensional function Ax 2 
with 1(9). For the above parameters, this yields in the form ^^j^: 



r 



V 2 . 
A mm 



a 

b 

M 



+0.0003, 
-0.0003, 



-0.0009 
-0.0009' 



-0.003,+0.009 
-0.003,-0.009' 



+0.02, +0.05 1f) 6 p_ v 
-0.02,-0.05 ' iu ot; v • 



(30) 



These are the ranges that we will use in the subsequent analysis. Note that in this model, 
they are already quite narrow, which is a manifestation of the fact that this model holds a 
lot of structure. 



3.3.3 Predictions for various LFV-experiments 

The most important types of LFV-experiments are rare lepton decays, — > ej-7, as well 
as conversions of a bound muon to an electron for some nucleus N, fiN — > eN. In a 
Ma-like model, the decisive quantities for both types of processes are [38J (ij = — ► 
e^/ej-ey-conversion): 



a 



ij 



2m 2 (H±) 



k=l 



[rrii + mj)I a 



Ml 
m 2 (iY±) 



+ M k I b 



Ml \ 
m 2 ( J ff ± ) ) 



(31) 
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Experiment 


Status 


Process 


BR-Limit / Sensitivity 


MEGA 


Past 


U — > e7 

r 1 


1.2 • 10 -11 


MEG 


Future 


a e7 

r i 


1.0 • 10" 13 


BELLE 


Past 




4.5 • 10" 8 


Babar 


Past 


t — »• e7 


1.1 • 10" 7 


MECO 


Cancelled 


/xAl — >■ eAl 


2.0 ■ 10" 17 


SINDRUM II 


Past 


/iTi — > eTi 


6.1 ■ 10" 13 


PRISM/PRIME 


Future 


/iTi — > eTi 


5.0 ■ 10" 19 


SINDRUM II 


Past 


liAu — > eAu 


7.0 ■ 10" 13 


SINDRUM II 


Past 


/iPb -> ePb 


4.6 ■ 10 -11 



Table 5: Limits on the branching ratios for several past and future LFV-experiments [9]. 



where 

Ia(t) 



16tt 2 



~2t 2 + 5t- 

. i2(t-i; 



t 2 Int 



2(t 



and hit) 



16tt 2 



2{t 



tint 



Using these, the branching ratios for the processes are given by 



Br(e ): 



ejj) 



3 

mj 



<7i 



8vr rfe,- 



-— and Bt(liN -> eN) 



D 2 



2 5 m 2 ^ capt (A0 



\o 



f ie | 



(32) 



(33) 



In the first formula, we have neglected the final state lepton mass. The quantities and 
a; capt (A r ), as well as a general expression for the second formula are given in Ref. [39] . 



3.3.4 Past and future LFV-experiments for Model 1 

We then use the parameter ranges from Eq. (1301) to make predictions with Eq. (133]) . The 
result is included in Fig. [TJ Furthermore, we have put in the limits/sensitivities of several 
past/future experiments, all listed in Tab. [5j A further discussion of the results will be 
given in the next section. 



4 Results 

We will now discuss how the general conflict between an extended scalar sector and flavour 
symmetries looks in our example models. Let us first start with Model 1. The numerical 
results can be seen in Fig. [TJ On the left panel, we present the la (black) and 3a (gray) 
predictions of Model 1 for the processes \i — > ej, r lij, and r — > ej, as well as different 
present and future bounds from several experiments, cf. Tab. [5j The right panel shows the 
same for ii-e conversion on the elements Al, Ti, Au, and Pb. 

Model 1 is the prime example that our logic works: As explained in Sec. 13.2.11 there 
are only 3 free parameters in the model. Still, it is able to fit the neutrino data well. 
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Figure 1: The numerical results of our analysis for Model 1. 



Actually, the only deviations from a perfect fit arise from the very accurate prediction 
of the mixing angles (e.g., the experimental best- fit value of #23 is not exactly maximal; 
cf. Eq. (1151) and Tab. H]). The obtained parameter ranges are, however, quite narrow, as 
can be seen from the example given in Sec. 13.3.21 This is exactly the point, where the 
experimental limits on LFV-processes get really powerful: Because of the stiffness in the 
model parameter space, the prediction of, e.g., the branching ratio /1 — > is so clear, 
that only a very narrow window is left for parameter variations. Accordingly, this model 
is actually already excluded by the past MEGA experiment (cf. Fig. [1]) for all four Higgs 
scenarios from Eq. ff27|) . We want to stress again, that these four scenarios belong to the 
few regions in parameter space that are indeed consistent with all the data and constraints 
mentioned in Sec. 13.3.11 The branching ratios for fi-e conversion are in general lower, and 
pass all current constraints. However, in this sector PRISM/PRIME will provide another 
future bound that will be able to exclude this model. 

The remaining questions is how far we can stretch this logic for models with less and less 
predictivity. As example for that case we can use Model 2, which has seven free parameters 
to fit the data (cf. Sec. 13.2.21) . This more than doubles the degrees of freedom in the fit. 
The numerical results for this model are given in Fig. [21 First of all, it may look odd 
that here, all la and 3a regions are somehow narrow, except for r — > fi'y. This is simply 
because all branching ratios are essentially functions of the product \ab\ (where a and b are 
model parameters), while the one for r — > fi'-f is given by the sum of three contributions, 
which are proportional to \b\ 2 , |c| 2 , and |<i| 2 , respectively. This numerical example nicely 
shows how more freedom blows up the regions which are predicted by a certain model. 
Turning this argumentation round, a certain limit on some observable is weaker the more 
free parameters there are that influence the observable in question. 

However, even this model with much less predictivity than the one before can be excluded 
for some scenarios: Scenario 5 has already been excluded by the MEGA-experiment and 
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Figure 2: The numerical results of our analysis for Model 2. 



scenario 7 can be tested by MEG. This shows the strength of our considerations: Even 
for a model that has a lot of freedom our logic still applies in suitable settings, which are 
here given by the scenarios 7 and 5. Actually, even the scenarios a and (3 are not that far 
below the future MEG-bound, and especially a hypothetical future experiment aiming at 
t — »• /17 might be very suitable to exclude this particular model. 



5 Conclusions 

In this paper we have studied the conflict arising in models with an extended scalar sector 
and discrete flavour symmetries when confronted with LFV-bounds. We have illustrated 
this using two examples based on the Ma-model, one with an A 4 and with a D 4 symmetry. 
Since the first model exhibits a relatively rigid structure (only three free parameters), it is 
already excluded for all four scenarios by existing bounds. Even though the second model 
has more than twice as many free parameters, it can still be strongly constrained and two 
of the scenarios can either be excluded or tested in the near future. 

We want to stress, however, that our considerations are not at all restricted to Ma-like 
models, but should apply to a much wider class of theories. Models with a lot of struc- 
ture (meaning few parameters) may easily be excluded by existing or future LFV-bounds 
although they have no problems without the flavour symmetry. Even models with many 
parameters can at least be strongly constrained, if not excluded as well. 
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Appendix A: Group Theory of [18 



The group is a group which describes even permutations of four objects. It has two 
generators, 5* and T, that fulfill the relations 



S 2 = (ST) 3 = T 3 = 1. 



(A-l) 



The group has four inequivalent irreducible representations, 1, V, 1", and 3, which trans- 
form under the generators, S and T as follows: 



1 

r 



s = i, t = 1, 

S = l, T = u 2 , 
5 = 1, T = u, 


















2 




3 : 


H 





lu 2 









-1 


• 






VO 









K 2 


2 





(A-2) 
(A-3) 

° \ 2 2 -1/ 

where u; = e l27r ^ 3 (which implies to 4 = to). 

The product rules for the singlets are the following: 

l'xl'^l", l'xl" = 1, l"xl" = l', 1x1 = 1, lxl' =1', lxl" = 1". (A-4) 

Consider now two triplets: 

a = (a u a 2 ,a 3 ) T , b = (&i, b 2 , b 3 ) T . 
The product of these two triplets can be decomposed as 

3x3 = l + l' + l" + 3. + 3 ( , I 



(A-5) 



(A-6) 



where 



1 = (afe) = ai&i + a 2 &3 + a 3 & 2 , 
1' = (a&)' = a 3 & 3 + ai& 2 + a2&i, 
1" = (ab)" = a 2 b 2 + a ± b 3 + ash, 



(A-7) 



and 



(a&) a = -(2ai&i - a 2 & 3 - a 3 & 2 , 2a 3 6 3 - ai& 2 - a 2 b x , 2a 2 b 2 - a x b 3 - a 3 b{f , 



3 a 



(a6) a = -(a 2 b 3 - a 3 b 2 , a x b 2 - a 2 h, a x b 3 - a 3 &i) . 



(A- 
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Appendix B: Group Theory of D 4 fl3l 140 



The group -D4 is a group which describes the symmetry of a square. It has two generators, 
A and B, that fulfill the relations 



A = B = 1 and ABA = B. 



(B-l) 



The irreducible representations consist of four singlets, Ii,l2>l3'l4' an( ^ one doublet 2, 
which transform under the generators, A and B as follows: 



±1 


A 


= 1, B = 


1, 




A 


= 1, B = 


-1, 


Is- 


A 


= -1, B 


= 1, 




A 


= -1, B 


= -1, 


A = 


(0 


*-*) • B 


-(! 



(B-2) 
(B-3) 



The product rules for lj are 



li x li = Il , Ii x lj = lj for i = 1, 4 , I 2 x 1 3 = I4 , 1 2 x 1 4 = 1 3 , and 
l3 x 14 = 12 (B-4) 



For Si ~ lj and (ai, a2) T ~ 2 we find 

^ 2 , ( 21 |~2 

-s 2 a 2 



sia 2 



s 3 ai 



and 



s 4 a 2 
— S4CL1 



2 . (B-5) 



For (ai, a 2 ) T , (61, fr 2 ) T ~ 2, the product 2x2 decomposes into the four singlets which read 
aifr 2 + a 2 6i~l_l , aife 2 — a 2 6x ~ 1_2 , ai&i + a 2 5 2 ~ 1_3 and ai&i — a 2 6 2 ~ 1_ 4 . (B-6) 
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